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A non-modal finite element theory for predicting the attenuation of an extended-reacting liner containing 
a porous facesheet and located in a no-flow duct is presented. The mathematical approach is to solve separate 
wave equations in the liner and duct airway and to couple these two solutions by invoking kinematic constraints 
at the facesheet that are consistent with a continuum theory of fluid motion. Given the liner intrinsic properties, 
a weak Galerkin finite element formulation with cubic polynomial basis functions is used as the basis for 
generating a discrete system of acoustic equations that are solved to obtain the coupled acoustic field. A state- 
of-the-art, asymmetric, parallel, sparse equation solver is implemented that allows tens of thousands of grid 
points to be analyzed. A grid refinement study is presented to show that the predicted attenuation converges. 
Excellent comparison of the numerically predicted attenuation to that of a mode theory (using a Haynes 25 
metal foam liner) is used to validate the computational approach. Simulations are also presented for fifteen 
porous plate, extended-reacting liners. The construction of some of the porous plate liners suggest that they 
should behave as resonant liners while the construction of others suggest that they should behave as broadband 
attenuators. In each case the finite element theory is observed to predict the proper attenuation trend. 
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= system matrix and diagonal matrix of density ratios, respectively 
= height of an element, width of an element, and sheet thickness, respectively 
= speed of sound 

= axial acoustic power and sound attenuation in dB, respectively 
= vectors of source effects and unknown node coefficients, respectively 
= height of duct airway and liner, respectively 
= unit imaginary number and wavenumber, respectively 
= location of leading and trailing edge of the liner, respectively 
= span wise mode order and 2D basis functions, respectively 
= number of elements in the duct airway and liner, respectively 
= number of grid lines in z and x directions, respectively, of duct airway 
= number of grid lines in z and x directions, respectively, of the liner 
= outward unit normal to a bounding surface and line, respectively 
= 3D and quasi-3D acoustic pressure field, respectively 
= bounding surface, volume, and bounding line, respectively 
= dimensional time and width of duct, respectively 
= 3D and quasi-3D acoustic particle velocity, respectively 
= transverse, spanwise, and axial Cartesian coordinate, respectively 
= real part of a complex quantity and vector dot product, respectively 
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v,V 
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UC 


= 3D and 2D gradient operator, respectively 
= mean density and unknown nodal coefficients, respectively 
= hole diameter and open area ratio, respectively 
= normalized acoustic impedance and impedance drop, respectively 


Subscripts: 

c,s,t,w 
dJA 
lx , /y, lz 
m,q 


= liner case, source plane, duct termination, and rigid wall, respectively 
= duct airway, facesheet, and liner quantity, respectively 
= effective density in the x, y, and z directions, respectively 
= spanwise mode order and basis function counter, respectivly 


Superscripts: 

* = complex conjugate 


I. Introduction 

S LIGHTLY more than four decades ago, the quest for more fuel efficient and quieter aircraft led to the replacement 
of the turbojet with the turbofan engine as the dominant propulsion system for wide-bodied commercial aircraft. 
Because fan-generated tones dominate the engine noise spectrum in turbofan engines, this provided the opportunity 
to exploit local-reacting acoustic liner technologies to suppress fan noise tones. These liners support air particle mo- 
tion normal to their surface and this motion depends only on the local acoustic field and normal incidence acoustic 
impedance. Local-reacting liners have been especially attractive in reducing turbofan engine noise because they pro- 
duce attenuation peaks in the vicinity of weakly damped resonances, and have therefore been ideally suited for the 
reduction of engine tones. 

The constraints imposed by a local-reacting liner are not sufficient for cases in which the motion at the surface of 
the liner depends on the joint action of the wave fields in both the liner and the duct airway. In these situations an 
extended-reacting liner model is essential. Extended-reacting liner structures normally consist of an airspace packed 
with foam or a fibrous material. These liners have generally been unacceptable for turbofan engine noise application 
because of the undesirable properties they exhibit within the extreme temperature and pressure environments of air- 
craft engines. Among the undesirable properties that render these liners non- flight worthy are fluid retention, debris 
contamination, wicking, lack of load bearing capacity, and degradation of sound absorbing properties. 1 

Over the last decade, two major changes to the turbofan engine design have generated renewed interest in extended- 
reacting liners for turbofan engine noise applications. First, the fan blade shapes were modified to incorporate lean 
and sweep, such that the acoustic energy in the blade passage frequency is rapidly attenuated. Second, the bypass 
ratio of turbofan engines continued to increase, with a resultant increase in thrust to drag ratio and fuel efficiency. 
These modifications caused the sound pressure level of the blade passage frequency to decrease relative to the other 
harmonics, while the broadband noise generated primarily by flow over the stators continued to increase. As a result, 
the most promising acoustic liners of the future are expected to attenuate broadband noise. The single-layer, local- 
reacting liner with its attenuation peaks in the vicinity of weakly damped resonances is clearly not suitable for this 
new role. Multi-layer, local-reacting liners provide acceptable attenuation over an extended frequency range, but suffer 
from additional cost and weight issues. 

Two recent developments that have the potential for broadband engine noise reduction are being investigated at the 
NASA Langley Research Center. The first is the use of metal foam liners. 2 The second development appears to offer 
broadband absorption while retaining the structural simplicity and mechanical advantages of locally-reacting liners 
and involves the implementation of perforations in the cell walls. 3-5 In order to take full advantage of these and other 
potential broadband absorption liner concepts for turbofan engine applications, appropriate extended-reacting models 
are required so that the liner can be designed for optimal noise suppression. 

While a number of methods are currently available for predicting the sound attenuation of extended-reacting liners 
in a duct, they have been generally directed to heating, ventilation, and air conditioning (HVAC) systems. 6-15 Their 
use for turbofan engine application is of questionable value for two primary reasons. First, they use fundamental mode 
analysis, which is not valid within a turbofan engine because the cross-sectional area distribution varies along the duct 
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axis. Second, these analyses were not designed for large diameter nacelles for which several hundred hard wall duct 
modes propagate without decay. 

The purpose of this paper is to describe ongoing research aimed at developing and validating a non-modal, finite 
element theory to predict the attenuation in flow ducts containing flightworthy, extended-reacting liners, such as those 
discussed in the previous paragraphs. For the sake of brevity, the decision was made to limit the scope of this paper to 
the no-flow effort. Results are presented showing how well the theory predicts the attenuation in ducts containing: (1) 
a Haynes 25 metal foam liner, and (2) an extended reacting liner with perforations in the cell walls. 

II. Problem Description 



Plane 


Figure 1. Three Dimensional Duct with extended-reacting liner and Cartesian coordinate system. 

This analysis considers acoustic wave motion in the no-flow, constant area duct depicted in Fig. 1. All transients 
are assumed to have propagated out of the domain or to have been dissipated away, so that only acoustic solutions 
of the form e l2n & are of interest. Inserted into the top wall of the duct is an acoustic liner with extended reacting 
material, whose purpose is to attenuate the sound. The liner is held within a rigid case for which the top wall, two 
sidewalls, front wall, and back wall are impervious to sound. A thin porous facesheet (not visible in Fig. 1) makes up 
the lower wall of the case. The acoustic properties of the facesheet are assumed described by a normalized impedance 
jump across the facesheet. Throughout this work, all impedances are assumed known, and are normalized by the 
characteristic impedance of the air in the duct airway. A right-handed Cartesian coordinate system, (x,y,z), with the 
z-axis pointing along the axis of the duct, the y-axis perpendicular to the duct sidewalls, and the v-axis perpendicular 
to the top and bottom wall is used. Consistent with most experimental apparatuses and aircraft engine nacelles, the 
source and termination planes are located in rigid wall sections of duct. As indicated in the figure, the source and 
duct termination planes are located at z = 0 and z = L, respectively. The source plane acoustic pressure is assumed 
known and the origin of the coordinate system is taken at the lower left corner of the source plane. The geometrical 
arrangement depicted in Fig. 1 is consistent with the NASA Langley Grazing Incidence Tube (GIT), Grazing Flow 
Impedance Tube (GFIT), and Curved Duct Test Rig (CDTR), where data will be measured to validate the current 
model. 

Generally speaking, the sound absorption of the material is a nonlinear function of the acoustic particle velocity 
vector, and methods of manufacturing and packing the material in the volume of space between the rigid case and 
facesheet may cause significant inhomogeneities in the material. Additionally, many materials are manufactured 
in such a manner that their elements align with a particular direction, creating acoustic anisotropy. Each of these 
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effects is accounted for in the current model. Further, the analysis makes no distinction as to whether the volume 
between the facesheet and rigid case is packed with a porous or fibrous material such as metal foam or is a cavity non- 
locally reacting liner without fill, such as a porous honeycomb liner. The analysis begins by describing the differential 
equations and boundary conditions that govern the acoustic field in the duct and in the liner. Kinematic constraints, 
consistent with a continuum theory of fluid motion, are then used to couple the sound field in the duct airway and the 
liner. 


III. Governing Equations, Boundary Conditions, and Kinematic Constraints 



Figure 2. Division of 3D duct into volumes and surfaces for computational analysis. 

Figure 2 is a schematic of the two regions of space to be analyzed. The composite structure consisting of the duct 
airway and liner is divided into volume, Vj, that is contained by the duct airway, and volume, V/, that is contained by 
the liner (V/ also contains the porous facesheet and liner casing). The duct volume, Vj, is bounded on the front by 
the source plane and in the rear by the termination plane. Part of the periphery of the duct is bounded by the rigid 
impervious left sidewall, right sidewall, lower wall, and upper wall. The remainder of the duct periphery contains the 
porous facesheet that is contained by the surface, Sf. The liner volume, V/, is bounded on the front, back, top, and 
sides by the impervious rigid walls of the case and on the bottom by the surface of the facesheet. With this division in 
mind, we begin our analysis by developing the equations for two separate acoustic boundary value problems. The first 
acoustic boundary value problem governs the solution in Vj, while the second governs the solution in V/. 

The acoustic boundary value problem in V c j has been discussed in some detail in Ref. 16 and is rewritten in compact 
vector format as 


=* =» 9 2nf 

V • V p + kjp = 0, k d = 

c, 7 

inVj 

0) 

II 

I^H 

on S s 

(2) 

-ik d p 
V /; • n = — — — 

=> V 

on S t 

(3) 

Vp» n =0 

on S w 

(4) 


where the termination impedance, is generally a function of the spatial coordinates x and y. The method used 
to determine the exit impedance is discussed in a later section of this paper. It is possible to replace Eq. (3) by a 
nonreflecting, non-local termination condition. The nonreflecting, non-local termination condition allows the sound 
field to exit the termination without reflecting and avoids the need to know the termination impedance. However, in 
the NASA Langley GIT, GFIT, and CDTR (where data to validate the current theory will be taken), experience to date 
shows that the duct termination is reflecting at some frequencies. For this reason, Eq. (3) is used. 
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The equation that governs linearized acoustic disturbances in the liner is 


v»[A,yvp+kjp = 0, [A/] 


Pd/Pix o 0 

o Pd/ Ply 0 

o o p d /p, z 



( 5 ) 


The acoustic waves in the airway [Eqs. (l)-(4)] are considered to be unaffected by viscous and heat conduction effects 
(i.e., the mean density, pj, and sound speed, cj, are real quantities). However, in the liner, viscous and heat conduction 
effects are important, and are modeled via the complex intrinsic properties (i.e., [A/] and c/) of the sound absorbing 
material. The physical justification for applying a Helmholtz equation similar to Eq. (5) to describe the acoustics 
within the extended-reacting liner was first proposed over sixty years ago by Scott 6 for isotropic and homogeneous 
acoustic materials. However, Eq. (5) accounts for both inhomogeneous acoustic materials and acoustic anisotropy via 
the effective complex density matrix, [A/], and effective complex sound speed, c/. Although these intrinsic properties 
are generally nonlinear functions of the acoustic particle velocity vector, they are assumed known in this analysis. 
Therefore, the governing differential equation [Eq. (5)] is considered linear. Several additional details concerning the 
derivation of Eqs. (5) may be found in the excellent paper by Peat and Rathi. 17 

Coupling the acoustic boundary value problem in Vj with that in V/ is necessary to obtain a unique solution for the 
acoustic pressure field. For this purpose, let p be decomposed as 


P = 


Pd in v d 
Pi in Vi 


( 6 ) 


Now a continuum theory requires that the normal component of the acoustic particle velocity remain continuous at the 
facesheet. When written in terms of the acoustic pressure field, this condition is 


[A /] V Pi • ~n = Vp d < 


on S f 


(7) 


The second condition to be satisfied by a continuum theory, is that the ratio of the pressure jump across the facesheet 
to the normal component of acoustic particle velocity equals the impedance jump across the facesheet. 15 When written 
in terms of the acoustic pressure field, this condition is 


-ik d { £ l _ Pd \ = AC on Sf 

v Pd • n 


( 8 ) 


Upon solving Eqs. (l)-(8) for the coupled acoustic pressure field, the acoustic power at an axial station, z, in the 
hard wall duct can be computed 18 


n // 1 r _> _> *i 

2 ^ Pd («<*•») dxdy, u d = - 


zXld 

i2nfp d 


(9) 


and the performance of the liner (i.e., the acoustic energy absorbed) is the decrease in decibels of the axial acoustic 
power from the source to the termination of the duct 18 


MB = 101og 10 




( 10 ) 


IV. Quasi-3D Formulation 

Consistent with the locally reacting theory, 16 the solution to the 3D problem is reduced to the solution of a series 
of 2D problems by assuming that the 3D solution may be represented as a series of quasi-3D hard wall duct modes 

?=£/>mCos(^) (11) 

This approach will be valid in a duct with rigid sidewalls, provided the liner does not scatter acoustic energy between 
the quasi-3D modes, p m . Implicit in this assumption is that the exit impedance, £*, for each quasi-3D mode is inde- 
pendent of the coordinate direction, y. It should be noted that the GIT, GFIT, and CDTR (with straight walls) satisfy 
the quasi-3D assumption. 
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I 



Figure 3. Slicing of 3D duct into surfaces and curves for quasi-3D analysis. 


The next task in the formulation is to take a cut through the duct in the (x, z) plane and define the appropriate 2D 
regions and bounding curves where the quasi-3D acoustic boundary value problems are solved. Figure 3 is a schematic 
of the two regions of space to be used in the quasi- 3D analysis. Sd and 5/ are the surfaces contained by the duct and 
liner, respectively, when the 3D duct is sliced parallel to the (x,z) plane. The 2D region is bounded on the front and 
back by the source boundary, denoted by T s and termination plane, denoted bt T t . The upper and lower boundary 
containing the rigid walls of the duct is denoted by T w , the rigid case is denoted by T c , and the portion of the upper 
wall containing the porous facesheet is denoted by Tf. 

The quasi-3D boundary value problem describing the acoustic field in Sd is now derived by substituting Eq. (11) 
into Eqs. (l)-(4) to obtain 


v.v 


Pd- 


kl- 


,ntn s ' 
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Pd 9 n = 


o 

II 

in S d 

(12) 

Pd = Ps 

on r. 

(13) 

-ikdPd 

on r t 

(14) 

• it = 0 

on r w 
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where the subscript m on pd has been dropped for the sake of brevity and p s , are understood to be functions of the 
span wise mode order, m, and spatial coordinate, x. Further, explicit expressions for the source pressure, p s , and exit 
impedance, £*, for selected examples are given in section VI of this paper. 

The corresponding boundary value problem in Si is 


V • [A;] 1 Pi + 



Pd 

P ly { w’ _ 


Pi = 0 , 


[A/] 


Pd/Plx 0 
0 Pd/Plz 
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(16) 


and the liner interface conditions are 


V pi »~ft = 0 
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(17) 
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( 19 ) 
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Figure 4. Finite element discretization of 2D duct and liner into x and z grid lines. 


V. Quasi-3D Finite Element Methodology 


The quasi-3D, finite-element method may be interpreted as an approximation of the acoustic field via an assem- 
blage of rectangular elements formed by the intersection of grid lines as illustrated in Fig. 4. It is assumed that, in the 
duct there are NZ grid lines in the axial and NX grid lines in the transverse direction. Similarly, the liner is discretized 
with nz grid lines in the axial and nx grid lines in the transverse direction as shown. A typical rectangular element, e , is 
shown in Fig 5. Each element consists of four local node numbers labeled 1, 2, 3 and 4, respectively, and is considered 
to have width a and height b. The finite elements need not be equally spaced, and to avoid violating the continuity 
and convergence requirements of the finite element discretization , 19 the constraint is imposed that no elements may 
straddle the facesheet. The objective of the method is to obtain the unknown quasi- 3D acoustic pressure at the nodes 
of each of the (NX — 1 ) (NZ — 1 ) elements in the duct and the (nx — 1 ) (nz — 1 ) elements in the liner. 

The finite element methodology uses C 1 continuous finite elements (i.e., the acoustic pressure and its derivatives are 
continuous at inter-element boundaries). Thus, within each finite element, p , is approximated by a linear combination 
of sixteen polynomials, Afi , A/ 2 , . . . Afi 6, which comprise a complete set of basis functions 


P = 


16 


q=\ 


( 20 ) 


3d d 2 

For the C 1 continuous elements, the coefficients <P q ’s are the values of p, and at the four nodes of the 

element, and the 2D basis functions, N q , are identical to those used in plate bending elements . 19 The sound source, 
dimensionless termination impedance, intrinsic properties of the material (i.e., [A /] and c/), and normalized facesheet 
impedance jump are represented in a series similar to Eq. (20), but the nodal values of these variables are assumed 
known. 


The weak Galerkin finite element formulation is implemented by substituting the trial solutions given by Eq. (14) 
into the differential equation governing each structure, and then rendering the resulting residual errors zero by requiring 
that they be orthogonal to the complete set of basis functions, N q 
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Figure 5. A typical 2D element and local node numbering system. 


where the sum in Eqs. (21) and (22) are to be performed using the rules of finite element assembly. The weak finite 
element Galerkin formulation then follows from the application of Green’s Theorem to the second derivative terms 
in Eqs. (21) and (22), followed by substitution of the hardwall or termination boundary conditions. For example, 
applying Green’s Theorem to the second derivative term in Eq. (21) gives 

J s (y.S?p d ) Nq dS = f s (V p d *V Nq^j dS - (y Pd .Jt)N q dT d (23) 

The weak form of Eq. (23) for an element along the exit boundary is obtained by substituting the exit plane boundary 
condition into Eq. (23) to obtain 


f s (y*?Pd)N q dS=J s (y Pd *?N q )dS + J r (24) 

However, the line integral in Eq. (23) vanishes along a rigid boundary, so that the weak form for elements that contain 
a rigid boundary is: 

j s (v.V Pd )N q dS=f s (? Pd .?N q )dS (25) 

The trial solution [see Eq. (20)] is chosen to satisfy the source and facesheet interface conditions by constraining the 
nodal degrees of freedom (<f>^) along the source and facesheet boundary. This leads to a large sparse, asymmetric 
matrix equation of the form 

[A]{d>} = {F} (26) 

A direct, asymmetrical, parallel sparse solver is used to solve for {d>} (i.e., the unknown nodal values of quasi-3D 
pressure and its derivatives at the nodes of the composite structure). All computations are performed on a parallel 
computer so that tens of thousands of degrees of freedom may be solved in minimal wall clock time. Finally, the 
right-hand side of Eq. (10) is integrated numerically, using the numerically computed sound field from Eq. (26) to 
calculate the sound attenuated in the duct. 
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VI. Results and Discussion 


The accuracy of the proposed finite element theory is demonstrated by simulating the attenuation of flightworthy 
extended-reacting liners located in a duct. The first illustrated example simulates the attenuation produced when a 
25.4 mm thick Haynes 25 metal foam liner is inserted into the GIT. The remaining illustrated examples simulate the 
attenuation produced by fifteen 50.8 mm thick porous plate liners inserted in the GFIT. Although efforts are underway 
to compare simulated attenuations from the current theory to measured values, these comparisons are outside the scope 
of this paper. 

A. Haynes 25 Metal Foam Liner Simulations 

A metal foam liner constructed from a cobalt alloy (Haynes 25) was selected for the first simulation. It was selected 
due to the fact that it has desirable properties when exposed to the extreme pressure and temperature environment of a 
turbofan engine. For the purposes of this investigation, it is assumed to be isotropic and homogeneous. Figure 6 shows 



Figure 6. Photograph of Haynes 25 metal foam samples. 

a photograph of three samples of Haynes 25 metal foam. Each sample has the following properties: 

1. Dimensions of 50.8 mmx50.8 mmx 10.8 mm 

2. Porosity of 3.15 pores per mm 

3. Density of 8% (i.e., 8% of the volume of a sample consists of the cobalt alloy) 

Two intrinsic properties (pi x and q) for the liner were measured in the Langley Research Center Normal Incidence 
Tube (NIT). A schematic of the NIT and supporting instrumentation used in the measurement is given in Fig. 7. 
Because Haynes 25 does not contain a facesheet, the impedance jump, A^, was set to zero in the numerical simulations. 
Recall that for an isotropic liner only one of the effective densities need be measured (i.e., pi x = p i y = p/ z ). 

The geometry of the airway used in this simulation is that of the NASA Langley Grazing Incidence Tube (GIT), 
for which L\ =203 mm, L 2 =610 mm, L = 813 mm, and H = W = 2h =51 mm. The GIT simulations are made at a 
static temperature of 297 K and static pressure of 101.3 kPa. The static density, pj, was computed from the ideal gas 
law. These values were used because they correspond to the conditions under which the liner intrinsic properties were 
measured in the NIT. It is noted that mode analysis can be applied to the GIT geometry, provided the hard wall sections 
upstream and downstream of the leading and trailing edges of the liner are removed. Therefore, the simulations are 
carried out in a section of the GIT between z = L\ and z = L 2 where single mode analysis can be applied directly to 
compare with the current finite element theory. 
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Figure 7. Schematic of Langley Normal Incidence Tube (NIT) with supporting instrumentation. 


The attenuation predicted from the finite theory is compared to that of a mode theory when the source is the least 
attenuated mode and the duct termination is anechoic. Such a comparison provides a valuable test case because the 
least attenuated mode is generally the most important in practical considerations. The modes for the extended-reacting 
liner were extracted using the approach of Ramakrishnan and Watson. 20 The source pressure, p s , corresponding to the 
least attenuated mode and normalized exit impedance, £*, corresponding to an anechoic termination were computed 
from the mode solution obtained from Ref. 20 and then used as inputs to the non-modal finite element theory. 



Figure 8. Comparison of attenuation spectra of least attenuated mode of Haynes 25 foam. 

The comparisons are presented in Fig. 8 for a 254 mm long (/ = 254 mm) sample of liner. Predictions are made 
from 0.4 to 3.0 kHz in increments of 25 Hz. These correspond to frequencies for which only the plane wave mode 
propagates without decay in the hard wall sections of duct. The non-modal finite element predictions were obtained 
on an evenly spaced 31 x 12 grid, an evenly spaced 16x7 grid, and an evenly spaced 7x4 grid. The non-modal finite 
element theory and mode theory are in excellent agreement on each of the three finite element grids. It is especially 
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surprising that on the coarsest grid (i.e., the 7x4 finite element grid), the discrepancy between the modal and non- 
modal theory is barely noticeable. Non-modal finite element results were run on an Altix work station using a single 
CPU, and a typical run on the finest grid (the 31x12 finite element grid) consumed less than 0.1 seconds of wall clock 
time. Although the non-modal finite element code is capable of extracting the solution at tens of thousands of grid 
points, this option was not required because of the small dimensions of the GIT. 

B. Porous Plate Liner Simulations 



Figure 9. Schematic of a porous plate liner. 

Twenty-seven porous plate liners were initially considered for this study. Each porous plate liner consisted of 
a series of porous partitions arrayed in the z and v dimensions as shown in Fig. 9. For convenience, each porous 
partition used within a liner was assumed to be identical and the liner is assumed isotropic and homogeneous in these 
two dimensions. Three geometric parameters were used to define these porous partitions; sheet thickness (x), hole 
diameter (d), and open area ratio (a). These parameters were assumed to take the following values: 

1. sheet thicknesses (x) of 0.016, 0.032, and 0.048 inches 

2. hole diameters (d) of 0.016, 0.032, and 0.048 inches 

3. open area ratios (a) of 0.045, 0.09, and 0.135 
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Given the sheet thickness, hole diameter, and open area ratio, an in-house computer code (ZKTL) was used to predict 
the intrinsic properties of each of the twenty-seven liners for frequencies ranging from 0.4 to 3.0 kHz (in 100 Hz 
increments). Only fifteen of the twenty-seven liners were observed to have distinct intrinsic properties. These distinct 
intrinsic properties were then used as inputs to the non-modal finite element code and the attenuation from each of the 
fifteen liners was predicted. 

The fifteen porous plate liners were numbered according to the value of porous partition resistance (e.g., liner #1 
contained the lowest porous partition resistance, liner #2 has the second lowest porous partition resistance, liner #3 
has the third lowest porous partition resistance, etc). Because the liners with the lower numbers (i.e., liner #1, liner #2, 
etc.) have the lowest porous partition resistances, they are expected to behave as resonant liners (that is, they should 
exhibit large attenuations at a few selected frequencies in their spectrum, but the liner attenuation is expected to drop 
sharply away from these selected frequencies). On the other hand, liners with the higher numbers (e.g., liner #13, liner 
#14, and liner #15) have the highest porous partition resistances, and are expected to behave more as broadband liners 
(that is, they are expected to exhibit more broadband attenuation and may therefore be preferred for the next generation 
of aircraft engines). The purpose of the simulations presented in this section is to verify that the finite element theory 
predicts these expected trends. 

The geometry of the duct is that of the NASA Langley Grazing Flow Impedance Tube (GFIT), for which L\ = 357 mm, 
L 2 =660 mm, L - 102 m, H =64 mm, h = 51 mm, and W = 51 mm. The GFIT simulations are made at a static temper- 
ature of 297 K and static pressure of 101.3 kPa. The effects of the hard wall section upstream and downstream of the 
liner are included in these simulations. Further, because the frequencies under consideration are below those where 
higher order modes may propagate in hardwall sections without significant decay, the source was set as a plane wave 
(p s = 1) and a nonreflecting plane wave termination was used at the exit plane (i.e., Cj ml). 

Attenuations for frequencies ranging from 0.5 to 3.0 kHz were predicted for each porous plate liner on four finite 
element grids, Gl, G2, G4, and G8. The G1 grid used a 31 x 4 evenly spaced grid in the duct airway and a 11 x 4 
evenly spaced grid in the liner, the G2 grid was twice as dense as the Gl grid, the G4 grid was twice as dense as the 
G2 grid, and the G8 grid was twice as dense as the G4 grid. Fig. 10 shows the results of the grid refinement study on 



Figure 10. Grid refinement study of the finite element method for liner # 1 on evenly spaced grids (Gl = 31x4 grid in the duct and 11x4 
grid in liner, G2=62x 8 grid in the duct and 22 x 8 grid in the liner, G4=124x 16 grid in the duct and 44 x 16 grid in the liner G8=258x 32 
grid in the duct and 88 x 32grid in the liner. 

liner # 1. Note that the predicted attenuation has converged on the G2, the G4, and the G8 grid. Similar convergence 
trends were observed for all fifteen liners although these results are not presented for the sake of brevity. 
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(a) liners #1 -liner #5 


(b) liners #6-liner #10 



(c) liners #1 1 -liner #15 


Figure 11. Predicted attenuation of extended reacting liners. 


Figure 11 compares the attenuations for porous plate liners #1 - # 15 using the densest grid (i.e., the G8 grid). 
As expected, the finite element theory predicts that the liners with the lower porous partition resistance (e.g., #1- 
#5) behave as resonant liners and the liners with the higher porous partition resistance (liners #10-# 15) behave as 
broadband attenuators. Because the next generation aircraft engines are expected to have significant energy in the 
broadband engine noise spectrum, the designs that produced the lower numbered liners (such as liner #1) may not be 
desired. On the other hand, liners with the higher numbers (such as liner #15) exhibit more broadband attenuation and 
may be preferred for the next generation of aircraft engines. 


VII. Concluding Remarks 

A non-modal, finite element theory, for predicting the attenuation in a no-flow duct with anisotropic, inhomoge- 
neous, extended-reacting liners containing a porous facesheet has been presented. Based upon the results of this study 
the following conclusions are drawn 

1. The grid refinement study shows that the predicted attenuation converges. 
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2. Attenuations predicted with this finite element theory are indistinguishable from those predicted using mode 
theory for a Haynes 25 metal foam liner (the intrinsic properties of the liner were measured in the Langley 
Research Center Normal Incidence Tube). 

3. Predicted attenuations for fifteen porous plate liners (some designed to behave as a resonant liners, and others 
designed to attenuate broadband noise) exhibit the expected trends. 

Results presented in this paper show that the theory predicts the appropriate attenuation trends and that predicted 
attenuation spectra are now ready to be compared with measured values. 
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